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ABSTRACT 

Context. The spin rate of stars evolves substantially during their lifetime, due to the evolution of their internal structure and to external 
torques arising from the interaction of stars with their environments and stellar winds. 

Aims. We investigate how the evolution of the stellar spin rate affects, and is affected by, planets in close orbits, via star-planet tidal 
interactions. 

Methods. We used a standard equilibrium tidal model to compute the orbital evolution of single planets orbiting both Sun-like stars 
and very low-mass stars (0.1 Mq). We tested two stellar spin evolution profiles, one with fast initial rotation (1.2 day rotation period) 
and one with slow initial rotation (8 day period). We tested the effect of varying the stellar and planetary dissipation and the planet's 
mass and initial orbital radius. 

Results. For Sun-like stars the different tidal evolution between initially rapidly and slowly rotating stars is only evident for extremely 
close-in gas giants orbiting highly dissipative stars. However, for very low mass stars the effect of initial rotation of the star on the 
planet's evolution is apparent for less inassive (IM^) planets and for typical dissipation values. We also find that planetary evolution 
can have significant effects on the stellar spin history. In particular, when a planet falls on the star it makes the star spin up. 
Conclusions. Tidal evolution allows to differentiate the early behaviors of extremely close-in planets orbiting either a rapidly rotating 
star or a slowly rotating star. The early spin-up of the star allows the close-in planets around fast rotators to survive the early evolution. 
For planets around M-dwarfs, surviving the early evolution means surviving on Gyr timescales whereas for Sun-like stars the spin- 
down brings about late mergers of Jupiter planets. In light of this study, we can say that differentiating between one spin evolution 
from another given the present position of planets can be very tricky. Unless we can observe some markers of former evolution it is 
nearly impossible to distinguish the two very different spin profiles, let alone intermediate spin profiles. Though some conclusions 
can still be drawn from statistical distributions of planets around fully convective M-dwarfs . However, if the tidal evolution brings 
about a merger late in its history it can also entail a noticeable acceleration of the star in late ages, so that it is possible to have old 
stars that spin rapidly. This raises the question of better constraining the age of stars. 

Key words. Stars: rotation - Planets and satellites: dynamical evolution and stability - Planet-star interactions - 



1. Introduction 

The spin rate is an important quantity for the evolution of a star 
and also for the evolution of any planets orbiting close-in. 

The parameter that govems the direction of tidal evolution 
for a planet orbiting a star (or a satellite orbiting a planet) is 
the initial semi-major axis with respect to the corotation radius, 
the orbital radius where the orbital period matches the central 
body's spin period. For a planet interior to the corotation radius 
the planet's mean motion is faster than the primary's rotation, 
so the tidal bulge raised by the planet on the primary lags behind 
the position of the planet. The planet feels a drag force that slows 
it down and causes its orbital radius to shrink, in some cases 
leading to an eventual merger with the primary. However, for a 
planet exterior to the corotation radius, the tidal bulge on the star 
is in advance with respect to the position of the planet and tidal 
forces push the planet outward. 

The rotational evolution of Sun-like stars can be described in 
3 main stages: the pre-main sequence (PMS) stage, the zero-age 
main sequence (ZAMS) approach and the main sequence (MS) 



relaxation (Bouvier, 2008). During the PMS stage the young 
stars are observed to have a range of spin periods, typically from 
a few to ~ 10 days, and there is evidence that a highly efficient 
braking mechanism is at work (Herbst et al., 2007). It is still not 
clear what mechanisms are responsible for the observed distri- 
bution and the angular momentum loss, but these may be due to 
the interaction between the star and surrounding accretion disk 
(e.g., Ghosh & Lamb, 1978; Shu et al., 1994; Matt & Pudritz, 
2005b; Matt et al., 2010), powerful stellar winds (Hartmann & 
MacGregor, 1982; Hartmann & Stauffer, 1989; Tout & Pringle, 
1992; Paatz & Camenzind, 1996; Matt & Pudritz, 2005a, 2008; 
Matt et al., 2012), or other processes. Toward the end of the 
PMS phase, the fastest stars in the observed distributions appear 
to spin up in a way consistent with angular momentum conser- 
vation, while the rotation rates of the slowest rotators does not 
appear to change significantly. Thus, near the ZAMS, the spin 
period distributions are the widest, typically ranging from a few 
hours to ~ 10 days (e.g., Bouvier et al., 1997; Bouvier, 2008). 
Once on the main sequence, the stellar structure evolves slowly 
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enough that the torque from ordinary stellar winds becomes im- 
portant. Thus, on gigayear timescales, the average spin rates de- 
crease (Skumanich, 1972), and the range of observed spin rates 
narrows. To bracket the range of observed stellar spin rates, we 
consider here two populations: initially fast rotators, whose evo- 
lution follows the upper envelope of the observed spin rate dis- 
tributions, and slow rotators that follow the lower envelope. 

During the PMS and approach to ZAMS, the observed spin 
period distributions of M-dwarfs is qualitatively similar to that of 
Sun-like stars. Observations of young clusters constrain the ro- 
tation period of low-mass stars younger than a few xlOO Myr 
(Stassun et al., 1999; Herbst et al., 2001; Irwin et al., 2008) 
but that approach fails for old clusters due to the faintness of 
old M-dwarfs. Nonetheless, old slowly-rotating M-dwarfs have 
been detected (Benedict et al., 1998; Kiraga & Stepien, 2007; 
Charbonneau et al., 2009). Contrary to Sim-like stars that are 
mostly radiative except for a small (in terms of mass) convec- 
tive region at the surface, very low mass stars (M* < 0.35Mq) 
are entirely convective (Chabrier & Baraffe, 1997). For Sun- 
like stars with a radiative core, the interface between the core 
and convective envelope is thought to be important for the mag- 
netic dynamo, whereas in fully convective low mass stars other 
mechanisms have to be invoked to explain their observed mag- 
netic activity (Reiners & Basri, 2007). For example, Chabrier 
& Kuker (2006) showed that mean field modeling can produce 
a dynamo, which creates large-scale nonaxisymmetric fields 
and Browning (2008) showed that three-dimensional nonlinear 
magnetohydrodynamic simulations of the interiors of fully con- 
vective M-dwarfs can also produce a large-scale dynamo. 

The coupling between stellar spin history and tidal evolu- 
tion has been studied by Zahn (1994) for close binaries and by 
Dobbs-Dixon et al. (2004) for short period planets. Individual 
systems where tidal interactions are thought to have played a 
role have also been the subject of various studies. Lin et al. 
(1996) proposed that the planet orbiting 51 Peg stopped its disk- 
induced inward migration because of its presence outside coro- 
tation before the disk dispersal. Individual systems of the OGLE 
survey have been studied by Patzold et al. (2004). Some studies 
give constraints for stellar dissipation, such as Carone & Patzold 
(2007) for OGLE-TR-56b and by Lanza et al. (2011) for the 
CoRoT-11 system. 

In this study, we try to have a more general and system- 
atic approach of the effect of the stellar spin evolution on the 
tidal evolution of close-in planets. To this end, we couple stellar 
evolutionary models (Baraffe et al., 1998), wind parametrization 
(Bouvier, 2008; Irwin et al., 201 1) and tidal evolution. We con- 
sider two limiting cases for the stellar spin evolution that corre- 
spond to a star whose initial rotation is either very fast or very 
slow. These different evolutionary paths can be seen in Figure 1 
of Bouvier (2008) for Sun-like stars, or in Figures 13 to 15 of 
Irwin et al. (2011) for M-dwarfs . The slowly rotating stars be- 
gin with a rotation period of 8 days and the fast rotating stars 
with a period of 1.2 days. Both fast and slow rotators evolve 
as explained in Bouvier et al. (1997), where the loss of angular 
momentum due to the stellar wind is quantified given different 
star-dependent parameters among which is the rotation rate of 
the star. The higher the spin of a star, the more active the star, 
the stronger the winds and the stronger the braking. 

Here we use a standard equilibrium model to study the tidal 
evolution of planets orbiting stars. Our paper is structured as 
follows. The tidal and star evolutionary models are briefly dis- 
cussed in Section 2. Some preliminary analysis based on order of 
magnitude study are made in Section 3 before giving the results 
of the tidal evolution of planets around the two types of stars 



considered here in Section 4. Finally, in Section 5 we discuss the 
difficulty of Unking the results of tidal evolution and observa- 
tions and the important effect of late mergers on the rotation rate 
of the stars. 



2. Model description 

We have developed a model to study the orbital evolution of 
planets around stars by solving the tidal equations for arbitrary 
eccentricity and also taking into account the observed spin evo- 
lution of stars. 



2.1. Tidal model 

The tidal model that we used is a re-derivation of the equilib- 
rium tide model of Hut (1981) as in Eggleton et al. (1998). We 
consider both the tide raised by the star on the planet and by the 
planet on the star. We use the constant time lag model (Leconte 
et al., 2010) and the internal dissipation constant cr that was cali- 
brated for giant exoplanets and their host stars by Hansen (2010). 

Taking into account both stellar tide and planetary tide, the 
secular tidal evolution of the semi-major axis a is given by 
(Hansen, 2010): 



-^ = -^\Nam-^Na2(e)\ 
a dt Tp^ n J 

-^[NaHe)-—Na2{e)], 



where the dissipation timescale is defined as 



(1) 



1 



9MpiMp + M,)Rl0o-, 



(2) 



and depends on the stellar mass M., its dissipation cr. and the 
planet mass Mp. Qp is the planet's rotation frequency, and n is 
the mean orbital angular frequency. The planet parameters are 
obtained by switching the p and * indices. Nal(e) and Na2{e) 
are eccentricity-dependent factors, which are valid even for very 
high eccentricity (Hut, 1981): 



Nal{e) 

Na2(e) 



1 + 31/2e^ + 255/8e'^ + 185/16e^ + 85/64e^ 

(1 - e2)15/2 

l-i-15/2e2 + 45/8e4-H5/16e^ 
(1 - e2)6 ■ 



The secular tidal equations can be extended to arbitrary 
obliquity, which has been done in Leconte et al. (2010). The 
equations for the eccentricity and planetary rotation rate can be 
found in Bolmont et al. (201 1), where tidal evolution was studied 
for a planet-brown dwarf system. 

2.2. Planets model 

We consider planets with a wide mass range, from 1 to 5 Mj 
(Mj = mass of Jupiter). 

For terrestrial planets we use the dissipation values based on 
Earth's dissipation value. Neron de Surgy & Laskar (1997) in- 
ferred the quantity ^2,6^^® = 213 s from the DE245 data for 
Earth. k2,p is the Earth's potential Love number of degree 2, 
which is a parameter depending on the moment of inertia of the 
body. It tells us how the body responds to compression (kj -3/2 
means that the body is an incompressible ideal fluid planet). ATp 
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is the constant time-lag. Hansen's CTp and the quantity k2,pATp 
are related through: 



k2,pATp 



2 ' 



(3) 



where Rp is the planetary radius and Q is the gravitational 
constant. 

The terrestrial planet's compositions are assumed to be an 
Earth-like mixture of rock and iron, following the mass-radius 
relation of Fortney et al. (2007). 

For 1 Mj, 2 Mj and 5 Mj planets, we use respectively the 
values of radius and mass of Jupiter, WASP-33b (Christian et al., 
2006; ColHer Cameron et al., 2010) and OGLE2-TR-L9 (Snellen 
et al., 2009). For the dissipation factor, we use Hansen (2010)'s 
estimates for gas giants : o-p = 2.006 x 10"^" g^^cm^^s"^ 

Because the planetary spin synchronization timescale is 
short compared to the other timescales considered here (Leconte 
et al., 2010; Heller et al., 2011), the planet rotation period is 
fixed to the pseudo-synchronization value at every calculation 
timestep. In this study, we chose not to treat the evolution of the 
obliquity of the bodies for simpUcity. 

2.3. Stellar evolution 

Sun-like stars and M-dwarfs are self-gravitating objects bom 
from a collapsing dense molecular cloud. In the beginning of 
their evolution they contract, and when the temperature at their 
core reaches the value of ~ 3 x 10* K, the PPI nuclear reaction 
starts and the stars enter the main sequence. The Sun is thought 
to have entered the MS around a few 10 Myr after its birth. After 
8 Gyr, the Sun will leave the MS to evolve towards a red giant. 

To compute the influence of dissipation into a star, one needs 
to know the internal structure, mainly the evolution of the ra- 
dius with time, the moment of inertia (through the gyration ra- 
dius, which is here considered constant) and the tidal dissipation 
factor cr». The first two quantities are provided by stars evolu- 
tionary models (Chabrier & Baraflte, 1997; Baraflre et al., 1998). 
However evolution models use unconstrained values for the ra- 
dius of the star at early ages (f < lO^yr), so some quantitative un- 
certainties can arise early in the stellar evolution (Baraffe et al., 
2002). However, here we consider evolution after fo = 5 Myr 
for Sun-like stars and after to = S Myr for M-dwarfs , so these 
uncertainties should remain negUgible. 

We assume that both Sun-like stars rotate as solid bodies, as 
in Bouvier et al. (1997), although more recent work has included 
the effect of internal differential rotation between the radiative 
core and the convective envelope (Bouvier, 2008). Given that 
low mass stars are believed to undergo the same kind of stel- 
lar wind braking as Sun-like stars (Irwin et al., 2011), we also 
consider sohd rotation (Morin et al., 2008). 

In this work, M, is held constant, and the eff'ect of mass loss 
(through processes Uke stellar winds) on the intemal structure of 
the star is considered negligible. 

2.3.1. Stellar dissipation 

The stellar dissipation factor is poorly constrained but Hansen 
(2010) gives estimates for Sun-like stars : cr, - 6.4 x 
10~55 g-icm'^^s'^cTj where cr» = 7.8 x 10 . Thus for Sun-like 
stars, the dissipation factor is cr, - 4.992 x 10"^* g"'cm"^s"'. 

Like a brown dwarf, a 0.1 Mq star is fully convective so we 
expect that the dissipation mechanisms within the O.IM© star 
should be closer to a brown dwarf than a Sun-like star. Thus, we 



use a dissipation factor of (XtjM = 2.006 x 10 *" g 'cm'^s"' for 
the O.lMg stars in our calculations (Hansen, 2010). No scaling of 
the M-dwarfs dissipation factor has been performed compared to 
the brown dwarfs value. In this work, we varied the stellar dissi- 
pation factor by a few orders of magnitude so that the real value 
for M-dwarfs is in the considered range. This dissipation is much 
larger than for a Sun-mass star, and as we will see in subsection 
3.2, this has implications for tidal evolution. Stars more massive 
than 0.35Mo that have radiative zones may be less dissipative 
and more similar to a Sun-like star. 

2.3.2. Rotational angular velocity 

Our calculations begin during the stellar PMS. The evolution of 
the observed spin distributions of PMS stars have often been pa- 
rameterized in terms of a "disk locking" scenario (e.g., Bouvier 
et al., 1997; Rebull et al., 2004, 2006; Edwards et al, 1993; Choi 
& Herbst, 1996), in which the spin period of the star is assumed 
to remain constant at a specified "initial" rate, for some specified 
amount of time (hypothesized to be associated with the dissipa- 
tion of the disk). For simplicity, we adopt this basic picture and 
start our calculations at the point of "disk dispersal," after which 
the stellar angular momentum evolution is computed according 
to physical equations. 

A primary goal of the present work is to determine how dif- 
ferent stellar spin histories influence the star-planet tidal interac- 
tion. To this end, we consider two diff'erent spin evolution tracks 
that approximately follow the fast and slow envelopes of the ob- 
served stellar spin distributions, and we do this for two diff'erent 
stellar masses, 0.1 and IM©. In order to describe these tracks in 
a physically self-consistent way, we adopt a simplified model for 
the stellar spin, adapted from Bouvier et al. ( 1 997) for solar mass 
stars and from Irwin et al. (201 1) for 0. IMg stars. 

For both 0.1 and IMq stars we add the effect of tides to the 
formula of Bouvier et al. (1997) for the loss of angular momen- 
tum due to the stellar winds (based on the formulae of Kawaler 
(1988) and MacGregor & Brenner (1991)). The expression for 
the angular momentum loss rate is: 



-1 



1 dJ _ 
Jdt ~ 7 
1 h 



1/2 



-1/2 



-I- 



Nol(e) ■ 



No2(e) 



(4) 
(5) 



where h is the orbital angular momentum, n is the mean orbital 
angular frequency, is the stellar dissipation timescale, and the 
functions Nol and No2 are defined as: 



NoHe) = 
No2{e) = 



l-Hl5/2e^-H45/8e'^-H5/16e'^ 

(1 - e2)13/2 

(1 - e2)5 • 



Here K, and Lo^at are parameters of the model from Bouvier 
et al. (1997). To reproduce the present rotation of the Sun, we 
use the value oiK = 1.6 x lO''' cgs and ojsm = 14 Oq. Bouvier 
et al. (1997) showed that for fast rotators (O* > Osadt a =\ and 
for slow rotators (D* < to sat), a -2>. 

Irwin et al. (201 1) proposed various parameters to reproduce 
the observational data for stars with masses 0.1 < M/Mq < 0.35, 
and the spin evolution they simulated was calculated for a star of 
mass 0.25Mo. They also noted that the M-dwarf fast and slow 
rotators could not be fit with a single value of the parameter 
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K. We consider here that the spin evolution of O.IM0 stars can 
be described with the same parameters as O.25M0 stars, and we 
note that this assumption can explain the differences between the 
curves we show below and the curves of Irwin et al. (201 1). 
In this work, we used the values: 



OJsat - 0.65 Qq, 

Kfas, = 2.03 X lO'*^ cgs, 
K,,,„ = 8.0 X 10^5 cgs, 



for initially fast rotators 
for initially slow rotators 



These values allows us to reproduce the values of the spin of 
the star at f = 5 Gyr in the two extreme trends seen in Fig. 14 of 
Irwin etal. (2011). 



Radius of a 0.1 and a 1 M^^^ star and rotation period profiles 

0.1000 E 




10 10 10 
Age of star (yrs) 

Fig. 1. Radius (top panel) and spin (bottom panel) evolution of a O.IMq 
and a IMg star. In the bottom panel, the full and dashed dotted blue lines 
represent respectively the evolution of the rotation period of an initially 
fast rotating star and an initially slow rotating star with no planet. As 
in the top panel, the blue curves correspond to the O.IM0 star and the 
black curves to the Solar-mass star. The blue diamonds in the bottom 
panel correspond to the values of the spin at 5 Gyr for the two extreme 
trends of Irwin et al. (2011). The vertical dashed lines represent t(, for 
the two stellar masses. 



where fo corresponds to the time of disk dispersal, and ftuies is 
given by : 



_ 1 dQ.| 

J tides — ^ J. „ 

£2* at \R^=cst,rg2.^=cst 

= ^[NoUe) - ^No2(e)]. 



(8) 



Here, = jj^ is the ratio of orbital angular momentum h to 
spin angular momentum and rgl^ is the square of the parameter 
rgt (which is the radius of gyration of Hut (1981)) which is de- 
fined as : / = Mt(rgtRt)^, where / is the moment of inertia of 
the star /„,,„,, is given by : 



fwind ~ 



\-l/2 



-1/2 



if < (jJsat 

if Q* > ojmt 



(9) 



The integration of the equations of Section 2.1 was per- 
formed using a fourth order Runge-Kutta integrator with an 
adaptive timestep routine (Press et al., 1992). The precision of 
the calculations was chosen such that the final semi-major axis 
of each integrated system was robust to numerical error at a level 
of at most one part in lO^*. 

Fig. 1 shows that after the time of disk dispersal Tdisk (long 
dashed vertical lines in Fig. 1), the star spins up due to contrac- 
tion. After a few xlO^ yrs, stellar winds start to efficiently spin 
down the star. 



3. Order of magnitude analysis 

3.1. Parameters space 

We present results for 0.1 and \Mq stars'. We investigate the 
evolution of planets with masses of IM^, lOM®, lMuranus(= 
14.5M(g), IMj, 2Mj and 5Mj. For the moment we assume zero 
initial eccentricity so that only the stellar tide governs the evolu- 
tion. In Subsection 4.3.2 we present the outcome of an extreme 
case with an initial eccentricity of 0.01. 

We consider that planets begin their tidal evolution at the 
time of disk dispersal Tdisk- The planet formation timescale is 
proportional to the orbital frequency (Safronov, 1969; Raymond 
et al., 2007) and is thus far shorter than the disk lifetime for 
close-in orbits, but planet-disk interactions probably overwhelm 
planet-star tidal interactions during this time. For Sun-like stars, 
this time is taken to be : Tdisk - 5 x 10'' yrs, which then cor- 
responds to our initial time fo, and for M-dwarfs , fo = Tdisk = 
8 X 10'' yrs. 

We tested values for the stellar dissipation o", between 1 and 
1000 times the mean values given by Hansen (2010). 



Fig. 1 shows 0.1 and IMq star radius evolution track as well 
as the different spin period profiles. For both stellar masses, the 
slow rotators have an initial peiiod of f ,0 = 8 days and fast 
rotators have an initial period of P^o = 1.2 days, as in Bouvier 
(2008). 

Finally, the expression for the stellar rotation is : 



n.(f) = Q,(fo) X 



r§2,(fo) (Rtito) 



rg2,{t) \ Kit) 



• exp ftideAt^ X exp 1^ /,v,„rfdf j , 



(6) 
(7) 



3.2. Order of magnitude and timescales 

The stellar dissipation timescale is given by equation 2. This 
timescale depends on the stellar mass M», radius and tidal 
dissipation factor cr,, and the semi-major axis a of the planet. 
As Rt shrinks with time, it is clear that for fixed a the stellar 
dissipation timescale increases with time, so the stellar tide be- 
comes weaker. Assuming Mp <K M,, the ratio of the dissipation 
timescale of Sun-like stars (Tq) to M-dwarfs (T^m) is: 



' Some simulations were performed with a stellar mass of O.8M0, but 
the results were very similar to what we obtained with the IMq star, so 
we do not discuss those results. 
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TdM 



Ram 



0"dM 
0"g ' 



(10) 



which gives the following values: 



10 TdM, 

TdM, 



1000 



at 8 Myr 
at 5 Gyr. 



We can see that comparing the magnitude of the tidal effect 
around a M-dwarf and a Sun-like star is equivalent to compare 
their radii and their dissipation factors. At all times. Sun-like 
stars have a bigger radius than M-dwarfs but have a smaller dis- 
sipation factor. 

At 8 Myr, the dissipation timescale of M-dwarfs is shorter 
than of Sun-like stars so the early tidal evolution will be 
more important around M-dwarfs than around Sun-like stars. 
However, at 5 Gyr, the dissipation timescale of M-dwarfs is 
much larger than of Sun-like stars so the late tidal evolution will 
be more important around Sun-like stars than around M-dwarfs. 



Evolution timescales 




0.007 



0.01 



0.02 0.03 
SMA (AU) 



0.04 0.05 0.06 0.08 



8 Myr and 5 Gyr the stellar dissipation timescale increases. At 
8 Myr, T, for a 0.1 Mq star is shorter than for a IMq star so the 
stellar tide will have a stronger effect on the planet orbiting an 
M-dwarf than a Sun-like star However, at 5 Gyr the stellar dissi- 
pation timescale of the M-dwarf is longer than 10 Gyr for planets 
at a > 0.007 AU. This means that for M-dwarfs the system will 
"freeze" after some time such that the interesting tidal evolution 
will only occur early in the evolution. For Sun-like stars, the in- 
crease of the stellar dissipation timescale is less pronounced so 
tides still matter for planets closer than ~ 0.018 AU at 5 Gyr. 

Concerning the evolution of t„ the trends are similar to those 
for r, but a few differences can be seen. In particular, the curves 
show a peak feature at a semi-major axis corresponding to the 
corotation radius. This is because Tg diverges when the semi- 
major axis is close to the corotation radius: if a planet forms 
precisely at the corotation radius of a non-evolving star, it will 
experience no tidal migration. The system is thus perfectly syn- 
chronized as the planet - which is already in synchronization - 
and the star always show each other the same sides. However, 
the corotation radius is an unstable equilibrium distance because 
inside the corotation radius - to the left of the peak - the planets 
migrate inwards and outside the corotation radius - to the right 
of the peak- the planets migrate outwards. For a system with an 
evolving star, the stellar tide always causes tidal migration. Fig. 
2 also shows that outward migration (taking place further away 
from the star) always occurs on longer timescales than inward 
migration (taking place closer in). 

Fig. 2 also shows that Jupiter-mass planets closer than 
0.02 AU around Sun-like stars with the mean dissipation value 
tidally migrate inward. Planets farther than 0.02 AU experience 
no noticeable semi-major axis change because t„ is so long. In 
section 4, we will show that Sun-like stars need a dissipation 
factor of 1000 x cr, in order to have noticeable tidally-induced 
changes within 10 Gyr. In Fig 2, this is equivalent to lowering 
the curves corresponding to IMq stars, effectively decreasing the 
different timescales so that interesting behavior can be observed 
within the 10 Gyr timescale we consider here. 

We can also consider the timescale of stellar spin evolution 
Tn, (again for zero eccentricity): 



Fig. 2. Stellar dissipation timescale T, and SMA-evolution timescale 
versus semi-major axis for stars of mass 0. IMg and IMq and a planet of 
Jupiter mass. The stellar dissipation timescale was calculated for mean 
dissipation factors and for two different star ages: 8 Myr and 5 Gyr. The 
difference between these two times is a difference in the radius of the 
star, which is smaller at 5 Gyr than at 8 Myr. 



Another important timescale is the semi-major axis evolution 
timescale (SMA-evolution timescale). Assuming zero orbital ec- 
centricity, it is given by: 



tAi- 



(11) 



In contrast with the stellar dissipation timescale r», t„ is af- 
fected by the spin of the star Tt is the limit of when — » 0. 

Figure 2 shows the dependance of the stellar dissipation 
timescale and the semi-major axis evolution timescale t„ on 
the semi-major axis a for a system with a Jupiter mass planet, 
for Mt = 0.1 M0 and IM©. As expected, the farther a planet the 
weaker the stellar tide. A Jupiter at a > 0.05 AU around a 8 Myr 
star of mass Q.IMq does not experience any noticeable semi- 
major axis evolution over 10 Gyr Fig. 2 also shows that between 



2n _ a 

7* \ n 



(12) 



where y, = is the ratio of orbital angular momentum h 
to the spin angular momentum and /» is the moment of inertia of 
the star. The ratio of the stellar spin synchronization timescale of 
Sun-like stars to M-dwarfs assuming Mp «c M* is: 



Q-dM rg2s 1 - OdM/» 
o-Q rgldM 1 - Oq/h ' 



(13) 



where rglj is the square of the parameter rg/ (which is the radius 
of gyration of Hut (1981)). 

This timescale depends on the moment of inertia of the star. 
If the star has a high inertia, the tidal forces need more time to 
bring the star into synchronization with the orbital frequency. 
Fig. 3 shows that the stellar spin evolution timescale is shorter 
for M-dwarfs than for Sun-like stars at 8 Myr, and it is longer for 
M-dwarfs than for Sun-like stars at 5 Gyr for close-in planets. 
Due to the strong dependence of tidal effect with respect to the 
orbital distance, the acceleration of the spin due to a planet inside 
the corotation radius is faster than its deceleration due to a planet 
outside the corotation radius. The spin-evolution timescale is 
short for close-in planets, when the planets fall towards the star 
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Fig. 3. Spin-evolution timescale versus SMA for stars of mass 0. 1 Mq 
and IMq and a Jupiter-mass planet. The spin-evolution timescale was 
calculated for mean dissipation factors and for stellar ages of 8 Myr and 
5Gyr. 



there will be an angular momentum transfer between the planets 
orbit and the spin of the star leading to a noticeable spin-up of 
the star. 

Stars generally have enough inertia that tidal interactions 
with a planet require longer than 10 Gyr to bring the star in 
synchronization. However, for a very dissipative M-dwarf and a 
Jupiter mass planet the stellar synchronization timescale is short 
enough to lead to a perfect synchronized system (as we will dis- 
cuss in Section 4.3.2). 

The timescale of planet spin evolution is usually much 
shorter than the stellar spin evolution timescale. That is why 
in this work we consider the planet is always in pseudo- 
synchronization -if the eccentricity is non zero- or in synchro- 
nization - if the eccentricity is zero. 

4. Results 

4.1. Orbital evolution of planets 

For Q.IMq stars, the effect of different spin profiles on the or- 
bital evolution of planets is apparent for Earth mass planets and 
mean dissipations values. For a IM^ planet beginning at an ini- 
tial semi-major axis of 9 x lO"-' AU, slow rotators systems will 
tend to make the planet fall on the star whereas fast rotators sys- 
tems will allow the planet to survive for 10 Gyr. 

Interesting effects start to occur for planets of mass higher 
than lOMg. Equation 2 shows that the stellar dissipation 
timescale also depends on the mass of the planet. The more mas- 
sive the planet, the shorter the dissipation timescale. For IM® 
planets, evolution timescales are too big to lead to significant 
changes in less than a few Gyrs. However for planets of mass 
higher than lOM®, the evolution timescales are compatible with 
visible changes over a few 10^ yrs. These planets begin to react 
as they cross the shrinking or expanding corotation radius. Fig. 
4 shows the results of simulations for lOM® planets orbiting a 
O.IM0 star, using mean dissipation values. A planet beginning 
at 9 X 10"^ AU falls on the star in a few hundreds of thousands 
years if the star is a slow rotator, but survives if the star is a fast 
rotator This latter planet experiences a small inward migration 



before its first crossing with the corotation radius and an outward 
migration after. For initial semi-major axes larger than 0.02 AU, 
the difference between the two spin profiles is negligible. One 
super Earth, GJ 1214 b, has been detected around a O.I6M0 star 
(Charbonneau et al., 2009), so this planet could be experiencing 
interesting tidal evolution. 
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Fig. 4. Tidal evolution of lOMe mass planets starting at different ini- 
tial semi-major axes around either a fast rotating or a slow rotating 
0. IMq star with mean dissipation factor. Top panel: evolution of the 
semi-major axis. The full colored lines correspond to fast rotating stars 
and the dash-dotted lines correspond to slowly rotating stars. The solid 
and dash-dotted red lines represent the evolution of the corotation ra- 
dius in both cases if there is no planet. The long black dashes represent 
the Roche limit. Middle panel: The corresponding stellar rotation evo- 
lution (the same line code is used). The black long dashes represent 
Tsat = '^^/(J^sai- Bottom panel : Equatorial velocity of the star vs time 
(the same line code is used). 



This effect is more dramatic if the stellar dissipation in- 
creased to 1 000 xcTj. The stellar tide contributes in pushing away 
the planets orbiting a fast rotating star very efficiently and we can 
find the same behavior as for massive brown dwarfs (Bolmont 
et al., 2011). For high dissipations and fast rotators, the planets 
are pushed faither away and for the inner most planets converge 
towards a given distance. If we considered an equally spaced dis- 
tribution of planets, the distribution would be more packed at the 
end of the evolution. 

As we saw in Subsection 3.2, the tidally-induced early semi- 
major axis evolution occurs on shorter timescale around a O.IM© 
star than around a IMq star So to be able to observe any discrep- 
ancies between spin profiles for Sun-like stars the stellar dissi- 
pation factor must be larger than the mean value. Hereafter, for 
Sun-like stars we use a dissipation factor of a thousand times 
Hansen (2010)'s mean value. 

Fig. 5 shows the tidal evolution of Jupiter-mass planets or- 
biting IMq stars. Differences can only be seen between fast and 
slow rotators for very close-in planets. For slow rotators, a planet 
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Fig. 6. Tidal evolution of a lOMg, mass planet starting at different initial 
semi-major axis around either a fast rotating or a slow rotating 0.1 Mq 
star. The full colored lines correspond to semi-major axis evolution. 
The full red line and the dashed red line respectively correspond to the 
corotation radius of a fast rotating M-dwarf and of a slow rotating M- 
dwarf with no planet. The black long dashes represent the Roche limit. 



Fig. 5. Tidal evolution of Jupiter mass planets of different initial orbital 
distance around either a fast rotating or a slow rotating IMg star with a 
large dissipation factor. Top panel: the colored lines represent the semi- 
major axis of the planets. Solid lines correspond to an initially fast ro- 
tating star and dashed dotted lines to an initially slow-rotating star. The 
red curves represent the corotation radius assuming no planet. Middle 
panel: evolution of the stellar rotation period. Bottom panel: Evolution 
of the equatorial velocity of the star. In middle and bottom panel, a red 
full circle represents the Sun's present rotation, and the red curves cor- 
respond to stars with no planets. 



beginning at 0.03 AU falls on the star in ~ 10'' yrs, whereas for 
the fast rotators it falls in < lO' yrs. 

The evolution of planets around fast rotators is interesting 
because planets beginning very close to the corotation radius 
experience an outward migration for a few xlO^ yrs. Around 
f = 50 Myrs, the star has spun down sufficiently due to stellar 
winds that these planets cross the outward-drifting corotation ra- 
dius. After the crossing, these planets migrate back inward and 
crash onto the star at aboutl Gyr. For planets beginning their 
evolution past 0.05 AU, the difference between the two spin pro- 
files is negligible. Planets beginning their evolution at a semi- 
major axis bigger than 0.08 AU survive the evolution on 10 Gyr 
timescales. 

4.2. Orbital evolution of planets compared with the age of the 
star 

The time on the x-axis in Figs. 4 and 5 starts at to, thus after 
5-8 Myr of stellar evolution, but the stellar evolution timescale 
may offer a more appropriate measure. Figure 6 shows the orbital 
evolution of a lOM® planet around a O.IMq star. We can see that 
the planets which do not survive crash on the star on a timescale 
much shorter than the age of the star. 

Planets orbiting O.IM© stars undergo similar evolution to 
planets orbiting brown dwarfs (see Bolmont et al., 2011). Most 
of the tidal evolution occurs at early times, when the radius of 
the stars/brown dwarfs is still large enough that the stellar/brown 
dwarf tide is strong enough to drive changes in the planets' semi- 



major axes. However, the radii of O.IM© stars and brown dwarfs 
decrease substantially in time such that after a certain interval the 
stellar/brown dwarf tide becomes weak and the system freezes 
(see explanation in Subsection 3.2). 

Thus, we expect that similar conclusions can be drawn for 
planets around O.IM© stars and brown dwarfs. A statistical dis- 
tribution of planets around fully convective M-dwarfs can pro- 
vide information about their dissipation factors. Indeed, as in 
Bolmont et al. (201 1), we can make inferences from how densely 
packed the orbital distribution of close-in planets is; a higher 
stellar dissipation leading to a more packed distribution of final 
semi-major axis, and increasing the shortest possible final semi- 
major axis (see Fig. 19 and 20 of Bolmont et al., 201 1). 

For Jupiter mass planets around IMq stars the situation is 
different. Planets which survived the early evolution due to the 
crossing of the shrinking corotation radius still fall on the star on 
Gyr timescale due to the spin-down of the star Figure 7 shows 
the evolution of Jupiter mass planets around a Sun-like star 

For fast rotators, either the planets fall very quickly com- 
pared to the stellar evolution or they fall in more than a few 
10^ yrs. Contrary to M-dwarfs , the tidal forces at late ages are 
still important for Sun-like stars (see Subsection 3.2). So the dif- 
ferences of evolution at early ages due to the difference of initial 
stellar spin period disappear in the end due to the long term evo- 
lution. No conclusions can be drawn as was the case for fully 
convective M-dwarfs . 

4.3. Planetary influence on spin evolution 

The spin evolution of a star with no planets is determined by its 
initial spin, its contraction rate and the efficiency of its stellar 
wind. However, if a planet is orbiting the star there are angular 
momentum exchanges between the orbit of the planet and the 
spin of the star In particular, if a planet is spiraling inward to- 
wards the star, it will have the effect of making the star spin up. 



7 



Emeline Bolmont et al.: Effect of the stellar spin history on the tidal evolution of close-in planets 



'+4-' P JL H 

^ ri^^' \- 




10 10° 
Age of star (yrs) 



10" 



Fig. 7. Tidal evolution of a IMj mass planet starting at different ini- 
tial semi-major axis around either a fast rotating (top panel) or a slow 
rotating (bottom panel) IMg star. The full colored lines correspond to 
semi-major axis evolution. The colored dashed lines correspond to the 
corotation radius and the red dashed lines correspond to the corotation 
radius of a star with no planet. The black long dashes represent the 
Roche limit. 
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Fig. 8. Tidal evolution of a Jupiter mass planet starting at different initial 
semi-major axis around either fast rotating or a slow rotating 0. 1 Mq star. 
Top panel: evolution of semi-major axis of a a planet around an initially 
fast rotating M-dwarf (bold solid lines) and around an initially slow 
rotating M-dwarf (dashed dotted lines). Bottom panel: Rotation period 
of initially fast rotators (solid lines) and slow rotators (dashed dotted 
lines). The black long dashes represent Tsm = In/aisat- The red curves 
represent the evolution of the spin of the star if there is no planet. 



4.3.1. Spinning-up and mergers 

Figure 8 shows the evolution of Jupiter-mass planets around ei- 
ther initially fast rotating or an initially slow rotating 0.1 Mq 
stars, using mean dissipation values. In Irwin et al. (2011) the 
time of the dispersal of the gas disk, Tdisk, varies with the na- 
ture of the star - fast rotator or slow rotator Here we considered 
that both fast and slow rotators decoupled from the disk at the 
same time fo = 8 Myrs. This way we can compare the effects 
of tides on the evolution of planets beginning in the same initial 
conditions - except of course the initial stellar spin. 

As above, we see in Fig. 8 that planets can survive in a 
wider range of initial semi-major axes around initially fast ro- 
tating stars. Compared with Fig. 4, the stellar spin evolution is 
more strongly modified by planetary migration. When the planet 
begins spiraling towards the star due to stellar tide, angular mo- 
mentum is transferred from the planet's orbit to the stellar spin. 
Thus, a falling planet can turn a slowly rotating star into a fast 
rotating star 

In this work we followed the prescription on Irwin et al. 
(2011) such that initially fast rotators always have Kf^st as 
wind parameter and initially slow rotators always have Kshw 
However, in the simulations from Fig. 8 some slow rotators 
actually become fast rotators early in their history. Given the 
prescription from Irwin et al. (2011), a slow- turned-fast rotator 
evolves back into a slow rotator due to the difference in wind 
parameterization. One can argue that a slow-turned-fast rotator 
should rather switch and spin down like a fast rotator instead. 



The late spin evolution of such a star would have little effect 
on any surviving planets because they must lie at larger orbital 
distances than the planet that perished and also because after a 
few xlO^ yrs, the radius of the star is small enough for the tidal 
effects to be negligible. 

The planet starting at 0.014 AU around a slowly rotating star 
(middle panel of Fig. 8) has an interesting evolution. The planet 
starts inside the corotation radius, so the stellar tide pulls it in- 
wards. Given the planet's large mass, it transfers a significant 
amount of angular momentum to the star and spins it up. The 
corotation radius thus shrinks until it catches up with the planet 
and inverts the tidal forces on the planet, which then experi- 
ences a slow outward migration. At 8 x 10^ yrs, the star has spun 
down due to the stellar winds sufficiently that the planet crosses 
back inside the corotation radius. For the rest of its evolution the 
planet migrated slowly inward. 

Stellar spin-up when a planet migrates in and falls on the 
star is also apparent for Sun-like stars. Figure 5 shows delayed 
peeks in the spin rate when the star experiences a merger^ with 
a Jupiter-mass planet. When a merger occurs, the star has an 

- Metzger et al. (2012) identified three different outcomes for the 
merger of a planet into a star, depending on at what orbital distance 
the Roche lobe of the planet becomes smaller than the actual size of the 
planet. They also discuss the observational signatures of such events. 
Some configurations lead to a luminosity peak which could last days 
or years. However, here we keep a simple description of merger events, 
and assume that a merger occurs when the planet reaches the Roche 
limit. Then, all the angular momentum is transferred from the planet's 
orbit into the spin of the star. 
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excess rotation which disappears on Gyr timescales due to stellar 
winds. For a system with a planet beginning at 0.06 AU, the 
merger occurs after 8 Gyrs, but just after the merger this old star 
spins five times faster than the Sun today, corresponding to a 
equatorial velocity of about 10 km/s. 

This effect is stronger for more massive planets. A 5 Mj 
planet beginning its evolution at 0.06 AU falls on the star in 
about 8 Gyr. This makes the parent star spin up to more than 
20 times faster than present Sun (see Fig. 9), with an equatorial 
velocity of about 60 km/s. This corresponds approximatively to 
the fast rotators' initial rotation period. 



^ Roche limit 
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Age of star - 1^ (yrs) 

Fig. 9. As Figure 5, but with 5Mj planets. We can see that the planet 
being more massive than in Fig. 5 the acceleration of the star is more 
visible. Some slow rotators even become fast rotators for a few 10 Myr 
because of the merger of the planet. When a merger occurs the star is 
spun up to a rotation rate which is of the order of 20 Q.q. 



Unlike close-in planets around O.IMq M-dwarfs , which sur- 
vived for Gyr timescales if they did not fall in first few Myrs, 
close-in planets orbiting Sun-like stars that survived the early 
evolution are still doomed to fall on the star This difference 
in behavior can be explained by the fact that M-dwarfs have 
much smaller radii than Sun-like stars at late ages and thus much 
weaker stellar tides. For a planet at 0.03 AU, the semi-major axis 
evolution timescale t„ for a 5 Gyr old M-dwarf is longer than 
10 Gyr (Fig. 2), but just 1 Gyr for a 5 Gyr old Sun-like star with 
a dissipation of 1000 x cr,. M dwarf stellar tides are simply not 
strong enough to make the planets fall at late times. 



4.3.2. Stellar synchronization 

We now study a system that displays relatively spectacular tidal 
evolution. The system is made up of a Jupiter-mass planet or- 
biting an initially fast rotating, highly dissipative (1000 x cr,) 
O.IMq star The planetary dissipation is set to zero to isolate the 
influence the stellar tide. 
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Fig. 10. Tidal evolution of a Jupiter mass planet around a fast rotating 
0. IMq star with a high dissipation factor. Top panel : Evolution of semi- 
major axis of Jupiter mass planets (full line) beginning at different initial 
semi-major axis and of the corresponding corotation distance (dashed 
line). The red dashed line corresponds to the corotation radius of a star 
with no planet. Middle panel : Evolution of the eccentricity of the plan- 
ets. Bottom panel : Evolution of the rotation period of the planets (full 
line) and of the corresponding star (dashed line). In both cases, the stel- 
lar spin initially increases to reach the spin period of the planet (we 
consider planets in pseudo-synchronization). Later in the evolution ei- 
ther the spin-up due to contraction or the spin-down due to stellar wind 
leads the spin of the star to differentiate from the spin of the planet. 



Tidal theory states that an equilibrium state is obtained for 
the two body problem when the orbit is circular, and when both 
orbit and spins are aligned and synchronized (Hut, 1980). In this 
equilibrium state, both planet and star thus have a spin which is 
equal to the mean orbital frequency, Qt - Q.p = n. The planet 
having a low inertia compared to the star, it will reach synchro- 
nization early in the evolution, typically in a few thousand years. 
However, usually the star because of its higher inertia never 
reaches synchronization. In most cases, it would require much 
more than 10 Gyr for the star to reach synchronization. However, 
in the example shown in Fig. 10, where the stellar dissipation has 
been significantly increased, the stellar tide is efficient enough to 
bring about synchronization of the stellar rotation with the spin 
of the planet in a few 10'* years to a few 10^ years. 

Fig. 10 shows a case. Let us focus on the case of the pur- 
ple curve which corresponds to a planet beginning at 0.014 AU. 
Initially the stellar rotation period is of 1 .2 days, which corre- 
sponds to an corotation distance of about 0.01 AU. In 10"^ years, 
the corotation radius moves outward and reached the planet's 
orbital distance (as can be seen in the bottom panel, the stellar 
spin period increases to match the planet's orbital period; re- 
call that we assume pseudo-synchronization for planets). After 
10^ years the stellar tide has made the eccentricity decrease to 
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zero. Indeed, because Q«/n < jj and because of the pseudo 
synchronization of the planet, this system would be in the region 
where e < of the phase diagram oj/nvs e presented in Leconte 
et al. (2010). Between 10^ years and 10^ years, the star and the 
planet remain in perfect synchronization with a circular orbit and 
spin synchronization. This is a equilibrium state for the system, 
or it would be if the star was not evolving in time. Indeed, the 
stellar radius shrinks which acts to spin-up the star while stellar 
winds act to spin-down the star, and these effects drive the sys- 
tem away from this equiUbrium state. After a few xlO^ years, 
contraction spins the star up, and then stellar winds make the 
star spin down until 1 Gyr The planet has crossed the corotation 
radius and the stellar tide causes it to fall slowly towards the star. 
The orbital angular momentum of the planet is transferred to the 
star and spins it slowly back up. 

The other (blue) planet in Fig. 10 shows a similar evolution 
but with two differences. First, the effect of the stellar tide is 
weaker because the planet is initially farther out (0.022 AU), so 
synchronization occurs later and for a shorter time. Second, the 
other physical phenomena influencing the spin evolution of the 
star kick in earlier. The acceleration phase is much more pro- 
nounced for this case than for the previous one. The second dif- 
ference is that initially, > jf j^e2(e) ' ^ phase 
space region where e > 0, so the stellar tide contributes in in- 
creasing the eccentricity. However, the spin of the star decreases 
faster than the mean orbital angular frequency n so around a few 
10* years, the system is back in the region where e < and the 
stellar tide contributes in decreasing the eccentricity. Between 
~ 6 Myr and ~ 2 Gyr, the planet is outside the corotation radius 
and Q,/n > so the stellar tide contributes in increasing 
the eccentricity once again. However, in the meantime the radius 
has shrunk to values such that the stellar tide is negligible and a 
small eccentricity is kept in the system. When the star starts to 
spin down due to the stellar wind, the both planet and star are 
once more in the phase space region where e < 0. 

This case represents an interesting case study although it is 
not likely that the stellar dissipation is strong enough to allow 
such evolution. 

4.3.3. Implication on stellar age determination 

As we have seen, tidal interactions between a star and a hot 
Jupiter can bring about significant acceleration of the spin of 
the star in some cases. This is the same main conclusion as Pont 
(2009), who explains that the observed excess rotation of stars 
with a transiting system can be due to the tidal interaction be- 
tween the planet and the star. We note that stars without planets 
may still bear the traces of their violent pasts due to planet-star 
mergers. We can see in Fig. 5 and 9 that while approaching the 
star, the spin of the star changes continuously to reach a maxi- 
mum at the merger. 

To determine the age of a star, one technique is "stellar gy- 
rochronology" (Barnes, 2003; Barnes et al., 2010; Meibom et al., 
2011; Epstein & Pinsonneault, 2012). It consists of measuring 
the spin of a star, and inferring the age of the star by assuming a 
spin history. Generally speaking, fast rotating stars are classified 
as young stars and slow rotating stars are classified as old stars 
that have been spun down due to stellar winds. 

In light of this study, we can see that this method might, in 
some cases, imply the wrong ages of rapidly rotating stars. A 
rapidly rotating star could be a young star but it could also be an 
old star which experienced a merger with a planet. For example, 
the merger of 5 Mj planet on an initially fast or slow rotating 



makes the star spin up almost to the initial rotation period of the 
group of fast rotators. Thus, old stars that are spun up would be 
mistaken for young objects. 

An old star that underwent a merger with a Jupiter size planet 
would reproduce several characteristics of young stars: they will 
be fast rotators, they will have infrared excess due to a hot dust 
disk and accretion signatures. A merger will also produce ex- 
treme UV signatures, soft X-rays and a peak of luminosity de- 
pending on the nature of the merger (Metzger et al., 2012). These 
stars would most hkely be mistaken for young stars. 

Gyrochronological ages of fast rotators should be checked 
against other techniques, in order to verify their youth. Another 
technique consists of estimating the age of the star by comparing 
its location on the Herzprung-Russel diagram with theoretical 
steUar evolution tracks. 



5. Conclusions 

We have shown that the stellar spin history affects the tidal evo- 
lution of close-in planets, albeit in a confined region of param- 
eter space. In cases where the stellar spin history does matter, 
the difference between two close-in planets - one orbiting an 
initially fast-rotating star and the other orbiting an initially slow- 
rotating star - comes from an early phase of outward tidal mi- 
gration for the planet around the fast rotator (caused by the star's 
closer-in corotation distance). This phase of outward migration 
does not occur for the planet orbiting the slow rotator. This out- 
ward migration and the decrease of the radius of the star weaken 
the later tidal evolution and effectively delay or sometimes even 
prevent the later in-spiralling of planets onto their stars. At later 
times both slow and fast rotators spin down due to stellar winds 
(Skumanich, 1972; Bouvier et al., 1997; Bouvier, 2008; Irwin 
et al., 2011), so the orbital history of close-in planets orbiting 
old stars depends on something that is not directly observable: 
the stellar spin evolution. 

For Sun-like stars the stellar spin history affects tidal evo- 
lution only in relatively extreme circumstances. In particular, 
the spin history has a strong effect if stars are very dissipative, 
with dissipation rates cr» of about 1000 times the fiducial value 
(Hansen, 2010). For strongly dissipative stars, there are differ- 
ences in tidal evolution for planets orbiting initially slowly- vs. 
rapidly rotating stars for very close-in (a < 0.05 AU), massive 
planets (Mp > Mj). 

In contrast, the stellar spin history plays a role in a much 
wider region of parameter space for O.IM© stars, mainly be- 
cause these stars are fully convective and so we think they 
are much more dissipative, with dissipation factors assumed to 
match those of brown dwarfs. For these stars the differences 
in tidal evolution for planets orbiting initially slowly- vs. fast- 
rotating stars are apparent for mean dissipation values, for plan- 
ets out to ~ 0.02 AU, and for planets with masses as small as 
IMe. 

Low mass stars and Sun-like stars have different dissipation 
factors and different radii, so the evolution timescales are differ- 
ent and evolve differently (see Figs. 2 and 3). The early tidal in- 
teraction is stronger for planets around very low mass stars, and 
the difference between fast rotator profile and slow rotator pro- 
file is apparent for planets of one Earth mass with mean dissipa- 
tion values. After some time, the system freezes in a given state 
because the radius of the star has shrunk too much for the tidal 
evolution to occur on less than 10 Gyr timescale. For Sun-like 
stars, the tidal evolution for mean dissipation factor occurs on 
very long timescales, which is why much more massive planets 
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and higher stellar dissipation rates are needed to produce stellar 
spin-driven differences in tidal evolution. 

The inward migration of a Jupiter orbiting inside the coro- 
tation radius of an initially slow rotating O-IM© star can lead to 
a significant spinning-up of the star. By the transfer of angular 
momentum from the planet's orbit, an initially slow rotating star 
can become a fast rotating star. This effect is more dramatic if 
the planet actually falls on the star. Irwin et al. (201 1) used dif- 
ferent wind parametrization for fast or slow rotators so they can 
infer from the present spin rate if the star was initially fast ro- 
tating or not. However we show here that it might not be that 
straight forward. If the star experienced a merger with a planet 
it can modify the rotation rate of the star and change the slow 
rotator into a fast rotator. 

Massive planets orbiting very low-mass stars with high dis- 
sipation rates (cr, x 1000) can create systems in perfect synchro- 
nization where the spin of the star is equal to the spin of the 
planet (Fig. 10). However, the equilibrium is not stable and the 
system departs from it as the star spins up due to contraction or 
spins down due to the stellar winds. This strongly alters the stel- 
lar spin profile because the star can be efficiently spun down by 
a planet initially located outside the corotation radius or spun up 
by a planet interior to corotation. 

Unfortunately, hot Jupiters around M-dwarfs are extremely 
rare due to the inefficiency of the planets formation processes 
around low mass stars (Laughlin et al., 2004; Ida & Lin, 2005; 
Kennedy & Kenyon, 2008). Only three hot Jupiters are known 
to exist around stars with masses less than Q.IMq (Pepe et al., 
2004; Hellier et al., 2011; Borucki et al., 2011; Johnson et al., 
2012) and none around stars with masses as small as O.IMq. Hot 
Jupiters around very low mass stars remain to be detected. 

A statistical distribution of planets around fully convec- 
tive M-dwarfs could constrain the tidal dissipation factor cr,. 
Specifically, from the location of the inner edge of the plane- 
tary semi-major distribution one can infer a inferior limit for the 
dissipation factor. The more distant the inner edge, the more dis- 
sipative the star. However, in order to draw these conclusions, 
one needs a good estimate of the stellar ages. 

For Sun-like stars such conclusions cannot be made because, 
if the dissipation rate is high enough to affect the orbital evolu- 
tion at early times, significant tidal evolution still takes place at 
late times, as well. Slow rotators and fast rotators have similar 
evolutions after a few 10^ yrs so the observation of a hot Jupiter 
orbiting a star of known age and known dissipation would not 
allow us to infer if the stellar spin history. Indeed, in both cases, 
different initial semi-major axis can lead to the same observed 
semi-major axis. One can imagine trying to infer a planet's or- 
bital evolution from the composition of its atmosphere to know if 
the planet came from a "cold" region (0.04 AU) or a "hot" region 
(0.03 AU). Unfortunately, this exercise is fraught with uncertain- 
ties in both the expected atmospheric composition of planets that 
form at different orbital distances and the tidal parameters (fi*,o> 
(T,, the age of observed stars). 

Nonetheless, we emphasize the planets crashing on the star 
at late ages can entail a significant spin-up of the star and cre- 
ate a population of old fast rotating stars. The spin-up of the star 
due to a merger has been pointed out in Levrard et al. (2009), 
where they found that planets falling on their host star due to 
tides never reach a tidal equilibrium. Jackson et al. (2009) also 
addressed the problem of tidally induced mergers and the effect 
of these mergers on the parent star. They also found that a con- 
siderable spin-up is to be expected and also a change in stellar 
composition. Planet-star mergers thus may confuse stellar age 
determinations. In general, fast rotators are thought to be young. 



although we have shown that a merger can lead to old, fast rotat- 
ing stars that would mimic many of the characteristics of young 
stars. An independent determination of the age of observed stars 
is therefore very important, especially for fast rotators. 

Acknowledgements. We are grateful to Andrew West for suggesting that we 
study low-mass stars rather than simply concentrating on the case of Sun-like 
stars. We thank Franck Selsis and Franck Hersant for their ideas and support. 
We also thank the CNRS's PNP program for funding, and the Conseil Regional 
d'Aquitaine for help purchasing the computers on which these calculations 
were performed. SPM acknowledges support by the ERC through grant 207430 
STARS2 (http://www.stars2.eu).. 



References 

Baraflfe, I., Chabrier, G., AUard, R, & Hauschildt, R H. 1998, A&A, 337, 403 
Baraflfe, I., Chabrier, G., AUard, R, & Hauschildt, P. H. 2002, A&A, 382, 563 
Barnes, R., Raymond, S. N., Greenberg, R., Jackson, B., & Kaib, N. A. 2010, 

ApJ, 709, L95 
Barnes, S. A. 2003, ApJ, 586, 464 

Benedict, G. R, McArthur, B., Nelan, E., et al. 1998, AJ, 116, 429 
Bolmont, E., Raymond, S. N., & Leconte, J. 2011, A&A, 535, A94 
Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 736, 19 
Bouvier, J. 2008, A&A, 489, L53 

Bouvier, J., Forestini, M., & Allain, S. 1997, A&A, 326, 1023 

Browning, M. K. 2008, ApJ, 676, 1262 

Carone, L. & Patzold, M. 2007, Planet. Space Sci., 55, 643 

Chabrier, G. & Baraffe, I. 1997, A&A, 327, 1039 

Chabrier, G. & Kuker, M. 2006, A&A, 446, 1027 

Charbonneau, D., Berta, Z. K., Irwin, J., et al. 2009, Nature, 462, 891 

Choi, R I. & Herbst, W. 1996, AJ, 1 11, 283 

Christian, D. J., PoUacco, D. L., Skillen, I., et al. 2006, MNRAS, 372, 1117 
Collier Cameron, A., Guenther, E., Smalley, B., et al. 2010, MNRAS, 407, 507 
Dobbs-Dixon, I., Lin, D. N. C, & Mardling, R. A. 2004, ApJ, 610, 464 
Edwards, S., Strom, S. E., Hartigan, P, et al. 1993, AJ, 106, 372 
Eggleton, P R, Kiseleva, L. G., & Hut, R 1998, ApJ, 499, 853 
Epstein, C. R. & Pinsonneault, M. H. 2012, ArXiv e-prints 
Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 
Ghosh, P & Lamb, R K. 1978, ApJ, 223, L83 
Hansen, B. M. S. 2010, ApJ, 723, 285 
Hartmann, L. & MacGregor, K. B. 1982, ApJ, 259, 180 
Hartmann, L. & Stauffer, J. R. 1989, AJ, 97, 873 
Heller, R., Leconte, J., & Barnes, R. 201 1, A&A, 528, A27+ 
Hellier, C, Anderson, D. R., Collier Cameron, A., et al. 201 1, A&A, 535, L7 
Herbst, W., Bailer- Jones, C. A. L., & Mundt, R. 2001, ApJ, 554, L197 
Herbst, W., Eisloffel, J., Mundt, R., & Scholz, A. 2007, Protostars and Planets V, 
297 

Hut, P 1980, A&A, 92, 167 

Hut, P 1981, A&A, 99, 126 

Ida, S. & Lin, D. N. C. 2005, ApJ, 626, 1045 

Irwin, J., Berta, Z. K., Burke, C. J., et al. 201 1, ApJ, 727, 56 

Irwin, J., Hodgkin, S., Aigrain, S., et al. 2008, MNRAS, 383, 1588 

Jackson, B., Barnes, R., & Greenberg, R. 2009, ApJ, 698, 1357 

Johnson, J. A., Gazak, J. Z., Apps, K., et al. 2012, AJ, 143, 111 

Kawaler, S. D. 1988, ApJ, 333, 236 

Kennedy, G. M. & Kenyon, S. J. 2008, ApJ, 673, 502 

Kiraga, M. & Stepien, K. 2007, Acta Astron., 57, 149 

Lanza, A. R, Damiani, C, & Gandolfi, D. 201 1, A&A, 529, A50 

Laughlin, G., Bodenheimer, P, & Adams, R C. 2004, ApJ, 612, L73 

Leconte, J., Chabrier, G., Baralfe, I., & Levrard, B. 2010, A&A, 516, A64+ 

Levrard, B., Winisdoerffer, C, & Chabrier, G. 2009, ApJ, 692, L9 

Lin, D. N. C, Bodenheimer, P, & Richardson, D. C. 1996, Nature, 380, 606 

MacGregor, K. B. & Brenner, M. 1991, ApJ, 376, 204 

Matt, S. & Pudritz, R. E. 2005a, ApJ, 632, L135 

Matt, S. & Pudritz, R. E. 2005b, MNRAS, 356, 167 

Matt, S. & Pudritz, R. E. 2008, ApJ, 681, 391 

Matt, S. P, Pinzon, G., de la Reza, R., & Greene, T. P 2010, ApJ, 714, 989 
Matt, S. P, Pinzon, G., Greene, T. P, & Pudritz, R. E. 2012, ApJ, 745, 101 
Meibom, S., Mathieu, R. D., Stassun, K. G., Liebesny, P., & Saar, S. H. 2011, 
ApL 733, 115 

Metzger, B. D., Giannios, D., & Spiegel, D. S. 2012, ArXiv e-prints 

Morin, J., Donati, J.-R, Petit, P, et al. 2008, MNRAS, 390, 567 

Neron de Surgy, O. & Laskar, J. 1997, A&A, 318, 975 

Paatz, G. & Camenzind, M. 1996, A&A, 308, 77 

Patzold, M., Carone, L., & Rauer, H. 2004, A&A, 427, 1075 

Pepe, R, Mayor, M., Queloz, D., et al. 2004, A&A, 423, 385 



11 



Emeline Bolmont et al.: Effect of the stellar spin history on the tidal evolution of close-in planets 

Pont, F. 2009, MNRAS, 396, 1789 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, 
Numerical recipes in FORTRAN. The art of scientific computing (Cambridge 

University Press) 
Raymond, S. N., Scale, J., & Meadows, V. S. 2007, ApJ, 669, 606 
Rebull, L. M., Stauffer, J. R., Megeath, S. T., Mora, J. L., & Hartmann, L. 2006, 

ApJ, 646, 297 

Rebull, L. M., Wolff, S. C, & Strom, S. E. 2004, AJ, 127, 1029 
Reiners, A. & Basri, G. 2007, ApJ, 656, 1121 

Safronov, V. S. 1969, Evoliutsiia doplanetnogo oblaka. (English transl.: 

Evolution of the Protoplanetary Cloud and Formation of Earth and the 

Planets, NASA Tech. Transl. F-677, Jerusalem: Israel Sci. Transl. 1972) 
Shu, F, Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781 
Skumanich, A. 1972, ApJ, 171, 565 

Snellen, I. A. G., Koppenhoefer, J., van der Burg, R. F. J., et al. 2009, A&A, 497, 
545 

Stassun, K. G., Mathieu, R. D., Mazeh, T, & Vrba, F J. 1999, AJ, 117, 2941 
Tout, C. A. & Pringle, J. E. 1992, MNRAS, 256, 269 
Zahn, J.-P. 1994, A&A, 288, 829 



12 



